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Abstract 
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1 Introduction 

Systems of singular differential equations have been a matter of increasing interest in theoretical 
physics and in some technical areas such as engineering of electric networks or control theory. 
The fundamental characteristic of these kinds of systems is that the existence and uniqueness 
of solutions are not assured. In particular, this situation arises in mechanics when dynamical 
systems described by singular Lagrangians are considered, and also when considering systems 
of PDE's associated with field theories described by singular Lagrangians (such as, for instance, 
electromagnetism) , as well as in some other applications related to optimal control theories. Fur- 
thermore, these systems do not have a nice Hamiltonian description, since not all the momenta 
are available and, in general, the equations have no solution everywhere. 

Bergmann and Dirac were pioneer in solving the problem for the Hamiltonian formalism of 
singular mechanical systems, by developing a constraint algorithm which gives, in the favourable 
cases, a final constraint submanifold where admissible solutions to the dynamics exist (in the 
sense that the dynamical evolution remains on this manifold) [SJ. Their main aim was to apply 
this procedure to field theories. Afterwards, a lot of work was done in order to geometrize this 
algorithm. The first important step was the work by Gotay et al |15j . and its application to the 
Lagrangian formalism [13} I14j. Other algorithms were given later, in order to find consistent 
solutions of the dynamical equations in the Lagrangian formalism of singular systems (including 
the problem of finding holonomic solutions) [3ll21|,H2]. and afterwards, new geometric algorithms 
were developed to be applied both in the Hamiltonian and the Lagrangian formalisms, as well as 
to other kinds of more general systems og singular differential equations [T71 [TU1 [2U1 E7J HD HI] ■ 

The Lagrangian and Hamiltonian descriptions of field theories is the natural extension of 
time-dependent mechanics. Therefore, in order to understand the constraint algorithm for field 
theories in a covariant formalism, the first step was to develop the algorithmic procedures for 
time-dependent systems. This work was provided mainly in EJ HU E7J EUJ ED ES ESI l38l [50] . 

There are several alternative models for describing geometrically first-order classical field the- 
ories. From a conceptual point of view, the simplest one is the k-symplectic formalism, which 
is the generalization to field theories of the standard symplectic formalism used as the geomet- 
ric framework for describing autonomous dynamical systems. In this sense, the /c-symplectic 
formalism is used to give a geometric description of certain kinds of field theories: in a local 
description, those theories whose Lagrangians or Hamiltonians depend on the fields and on 
the partial derivatives of the fields, or the corresponding moments, but not on the space-time 
coordinates [3D]. The foundations of the /c-symplectic formalism are the /c-symplectic mani- 
folds [Tj[2lE2]. Historically, it is based on the so-called polysymplectic formalism developed 
by Giinther [18] . who introduced the concept of polysymplectic manifold. Then, /c-symplectic 
manifolds are polysymplectic manifolds which have Darboux-type coordinates [32]. A natural 
extension of this formalism is the k-cosymplectic formalism, where /c-cosymplectic manifolds are 
used to describe geometrically field theories involving space-time coordinates or analogous ones, 
on the Lagrangian or the Hamiltonian [331 El] . This is the generalization to field theories of the 
cosymplectic formalism geometrically describing non-autonomous mechanical systems. One of 
the advantages of these formalisms is that one only needs the tangent and cotangent bundle of 
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a manifold to develop them. 

It is worth noting that G. Sardanashvily et al [121147] developed a polysymplectic formalism 
for classical field theories which differs from the one proposed by Giinther. (See also [22] for 
more details on the polysymplectic formalism.) In addition, we must remark that the soldering 
form on the linear frames bundles is a polysymplectic form, and its study and applications to 
field theory constitute the n-symplectic formalism developped by L.K. Norris [391 143] . 

Working within the framework of the A:-symplectic description for these theories, we present 
in this paper a geometric algorithm for finding the maximal submanifold where there are con- 
sistent solutions to the field equations of singular theories. This algorithm is a generalization 
of the presymplectic constraint algorithm for presymplectic dynamical systems [13], and gives 
an intrinsic description of all the constraint submanifolds. The problem is stated in a generic 
way for fc-presymplectic Hamiltonian systems, in order to give a solution to both Lagrangian 
and Hamiltonian field theories, as well as other possible kinds of systems of partial differential 
equations. In this framework, the solutions to these equations are given geometrically by inte- 
grable /c-vector fields in the manifold where the equations are stated. In this way, a constraint 
algorithm can be developed giving a sequence of submanifolds which, in the best case, ends 
in some final constraint submanifold where field equations have consistent solutions (fc-vector 
fields), although not necessarily integrable. The general problem of integrability is not addressed 
in this paper, only discussed in the examples. Finally, Lagrangian and Hamiltonian field theo- 
ries are particular cases where the above results are applied straightforwardly, although in the 
Lagrangian case the problem of finding holonomic solutions must be also analized. In addition, 
the unified Lagrangian-Hamiltonian formalism of Skinner-Rusk [32] , which was adapted recently 
for fc-symplectic field theories [45j, constitutes a framework where this algorithm is applied in a 
very natural way. A description of constraint algorithms for other geometrical models of field 
theories (multisymplectic) was made in pjjl [Til [23 [29] . 

The paper is organized as follows. In section 2 fe-symplectic structures are reviewed, as well 
as the corresponding Hamiltonian systems. Section 3 is devoted to fc-presymplectic Hamiltonian 
systems and the presymplectic constraint algorithm. In section 4 the particular case of field 
theories described by a Lagrangian function is considered, either in Lagrangian or in Hamiltonian 
formalism. Finally, the application to the Skinner-Rusk formalism and two examples are studied 
in section 5. 

Manifolds and maps are assumed to be smooth. Sum over crossed repeated indices is under- 
stood. 

2 £>symplectic Hamiltonian systems 

2.1 fc-symplectic manifolds. The bundle of A; 1 -covelocities 

Definition 1. A /c-symplectic structure on a differentiable manifold M of dimension N = 
n + kn is a family (w 1 , . . . ,uj k ; V), where each uj a is a closed 2- form, and V is an integrable 



X. Gratia, R. Martin, N. Roman-Roy — Constraint algorithm for k-presymplectic systems 4 

nk- dimensional tangent distribution on M such that 

(i) u A \vxV = for each A, 
(if) n^ =1 Kerw A = {0}. 

Then (M, u> A ,V) is called a £;-symplectic manifold. 

Theorem 1. fl\ \32$. Let (uj a ,V) be a k-symplectic structure on M. For every point of M there 
exists a neighbourhood U and local coordinates (q l ,p A ) (\ < i < n, 1 < A < k) such that, on U, 

uj A = dq i Adpf, y=/A... 5 JL 

\m m / i=x,...,n 

These are called Darboux or canonical coordinates of the fc-symplectic manifold. 

The canonical model of a /c-symplectic manifold is (T\,)*Q = T*Q(B k . ©T*Q, the bundle of 
k 1 -covelocities of an n-dimensional differentiable manifold Q, which has the natural projections 

7T A :(TlyQ - T*Q ; {T\fQ -» Q 

(q; a\, . . . , a%) ' * (q;a A ) ; (q;a], . . . ,a k ) h-> q 

(Tj,)*Q is endowed with the canonical forms 

e A = (7r A ye, u A = (tt a )*uj = -(ir A )*de = -de A , 

where 8 and to are the Liouville 1-form and the canonical symplectic form on T*Q. 

If (q l ) (1 < i < n) are local coordinates on U C Q, the induced coordinates (q l ,p A ) (1 < 
4 



A < k) on (7Tg) (U) are given by 



a],..., a k ) = q\q), p A {q; a\,..., a k ) = af 

Then we have 

A = p A dq\ u) A = dq { A dp A . 



d_ 

dq i 



Thus, the triple ((T\)*Q, uj a , V), where V = KerTvr^, is a /c-symplectic manifold, and the 

-I) 



natural coordinates in (T\,)*Q are Darboux coordinates. 



2.2 k- vector fields and integral sections 

Let T\M = TM© k . ©TM be the bundle of k 1 -velocities of a differentiable manifold M. It is 
endowed with the natural projections 

r A :T\M - TM ; r^: T\M - M 
(q,v lq ,...,v kq ) ' ^ (gjt'Ag) ; (g,uig,...,?;fc g ) i-" <7 

Definition 2. ^4 /c-vector field on a manifold M is a section X: M — > T^M of the projec- 
tion tIj. 



Therefore, giving a k- vector field X is the same as giving k vector fields X±, . . . ,X k on M, 
obtained as Xa = o X. We denote X = (Xi, . . . , X k ). 
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Remark The term fc-vector field on M is more often applied to the sections of the bundle 
A k TM — > M, that is, contravariant skew-symmetric tensor fields of order k. The A;- vector fields 
X = (Xi, . . . ,Xk) used here lead to a particular class of such tensor fields, the decomposable 
ones, X\ A ... A Xf~, which can be associated with distributions on M. 

Definition 3. An integral section of the k-vector field X = (X\, . . . , X^) is a map (ft: J — > M , 
defined on an open set J Cl' , such that 

d 

where t = (t 1 , . . . ,t k ) denote the canonical coordinates o/K fc . 
Equivalently, an integral section satisfies the equation 

cpW = Xo0, 

where (f>^ : J —> T^M is the first prolongation of <fi to T^M defined by 

* m W-(«t).T#(^|,),....T#(^| i )) . 

A k-vector field X is integrable if every point of M belongs to the image of an integral section 
o/X. 

d 

In coordinates, write Xa = X\—^. The <b is an integral section of X if, and only if, the 

ox 1 

following system of partial differential equations holds: 




Proposition 1. A k-vector field X = (X\,... ,Xj^) is integrable if, and only if, [Xa,Xb] = 
for each A,B. 

This is the geometric expression of the integrability condition of the preceding differential 
equation (see, for instance, or |SJ). 



2.3 Hamiltonian systems 

Definition 4. Let (M, uj a ,V) be a k-symplectic manifold, and a S f2 1 (M) a closed form. 
(M, oj a ,V, a) is said to be a fc-symplectic Hamiltonian system. 

As a is closed, for every point of M there exists a neighbourhood U C M and a function 
H € C°°(U) such that a = dH on U. This function is called a local Hamiltonian function. If a 
is exact, then H G C°°(M) is called a (global) Hamiltonian function. These functions are unique 
up to a constant on each connected component of M. From now on, we will write a = dH. 

The Hamilton-de Donder-Weyl (HDW) equation for a map ip: J — > M (J C is 

HS$)w A = dHoift . (2) 
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In canonical coordinates this reads 

dp _ dH dtpf dH 
W~dpJ' ~dt A ~~~W' 

where ip = (ifi 1 , ipf~). Recall that, according to our conventions, a sum is understood 
whenever the index A appears twice in upper and lower position. 

In order to give an alternative geometrical interpretation of these equations, we introduce the 
set 3Zjj(M) of those k- vector fields X = (X\, . . . , X^) on M which are solutions of the geometric 
field equation 

i{X A )u A = dH. (3) 

For fc-symplectic Hamiltonian systems, solutions of equation ([3|) always exist (this is a con- 
sequence of the lemma and the theorem in the next section). They are neither unique, nor 
necessarily integrable. 

d d 

In canonical coordinates of M, writing Xa = (XaY — — + (Xa)F „ R , equation (T3l) reads 

oq l dpf 

9H (V \A 9H (v \i 

W = ~ i A)l ' dpJ = {XA) - 

This geometric field equation for X is an alternative formulation of the HDW equation in 
the following sense: 

Proposition 2. Let X = (X\, . . . , Xk) be an integrable k-vector field in M. Every integral 
section ip: J — > M o/X satisfies the HDW equation (0|) if, and only if, X G X^(M). 

Note however that equations ^ and (J3l) cannot, in general, be considered as fully equivalent: 
a solution to the HDW equations may not be an integral section of some integrable ^-vector 
field on M. Solutions if) that are integral sections of some X G X k H {M) will be called admissible, 
and we will restrict our attention to them. 



3 £>presymplectic system and constraint algorithms 

3.1 /c-presymplectic Hamiltonian systems 

To consider singular field theories we have to drop some assumptions in the definition of a k- 
symplectic structure. So, a family (a; 1 , . . . ,uj k ) ok k closed 2- forms on a smooth manifold M 
will be called a k-presymplectic structure; accordingly, (M, uj a ) will be called a k-presymplectic 
manifold. 

The simplest example of a A:-presymplectic manifold is provided by any submanifold of a 
fc-symplectic manifold: the pull-back of the k 2- forms by the inclusion map yields k 2-forms on 
the submanifold. 

In some particular A;-presymplectic manifolds one can find Poisson-like coordinates, but it is 
an open question to characterize the necessary and sufficient conditions for these coordinates to 
exist. 
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Given a closed 1-form a G i? 1 (M), (M, to A ,a) is said to be a k-presymplectic Hamiltonian 
system. As above, we will write a = dH (locally or globally). 

Then we can also consider the Hamilton-de Donder-Weyl equation i(tj;^)uj A = dHoip, and 
also the set X k H (M) of fc-vector fields X that satisfy the geometric field equation i(Xa) oj a = dH. 

For fc-presymplectic systems the existence of solutions of this equation is not assured every- 
where on M. We will analyze the existence of solutions on a certain submanifold of M. 

3.2 Statement of the problem 

The problem we wish to solve arises from the Lagrangian and Hamiltonian fc-presymplectic 
formalisms in field theories, although other kinds of systems could also be stated in this way. 

Statement. Let (M,u> A ,dH) be a k-presymplectic Hamiltonian system. We want to find a 
submanifold C of M and integrable k-vector fields X = [X\, . . . , A&) G X k (M) such that 

i{X A )u A » dH (4) 

c 

( this means equality on the points of C) and 

X is tangent to C 
(this means that X\, . . . , X^ are tangent to C). 

As stated in the introduction, we will focus on the consistency of the equation and will not 
address the integrability condition in generality. 

Given a submanifold C of M, with natural embedding jc - C ^ M, let T k jc'. T\C — > 
T\M be the natural extension of jc to the /c-tangent bundles, and denote its image as T\C = 

We can define the map 

(p;v pi ,...,v Pk ) ' ^ (p,i(v Ap )u A ) 
and denote by (TC)j^ the annihilator of the image of T\C by bj,; that is, 

(TC)i = h(TlC)}° = {u p G TM | V(v pl , • • • , v p k) G T^C, i(v pA )u£ , u p ) = 0}. 

We call (TC)^ the k-presymplectic orthogonal complement of T\C in T\M. 
In particular, for C = M we have: 

fc 

Lemma 1. (TM)J = {(p,u p ) G TM | u p G Q Kerwp 4 } 

A=l 

(Proof) For every p G M and (u p i, . . . , f p ^) G (T£) p M, if u p G f] A =i Kerwp 4 , then we have 
(i(v Ap )oj A )(up) = -i(v A p)i{up)uj A = 0, and therefore u p G (T p M)^. 
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Conversely, if u p G (T p M)^, then (i(v Ap )^p)(u p ) = for every (v pl , . . . , v pk ) G {T\) p M. 
Then taking any (v p ±, 0, . . . , 0) with u p i / we conclude that u p G Ker w^; and analogously for 
the others. ■ 

The main result is the following: 
Theorem 2. Let C be a submanifold of M . The following conditions are equivalent: 

• there exists a k-vector field X = (X\, . . . ,X k ) G X k (M), tangent to C, such that equation 
l^p holds 

i(Y p )(dH) p = for every p £ C , Y p e (T P C)^. (6) 

(Proof) (=^) If there exists a fc-vector field X = (X\, . . . , X k ) G X k (M), tangent to C such 
that equation (JI|) holds, then, for every p£C and Y p G (T p C)^, 

= [i{X Ap )^]{Y p ) = i(Y p ) i{X Ap )^ = i(Y p ){&H) p . 

If ([6D holds, then 

(dfl) p G [(TC)tf = [h(TlC)] }° = h(T\C), 
and hence there exists (X p i, . . . ,X pk ) G (Tj,) p C such that (JH) holds. ■ 

3.3 /c-presymplectic constraint algorithm 

The application of the above result leads to an algorithmic procedure which gives a sequence of 
subsets . . . C Cj C . . . Ci C C\ C M. We will assume that: 

Assumption. Every subset Cj of this sequence is a regular submanifold of M . 

These submanifolds are sequentially obtained from the analysis of the consistency of a linear 
equation, namely eq. (|1J) at each point: 

i(X Ap )uj^ = (dH) p . 

First, C\ > M is the submanifold of M where this equation is consistent: 

d = {p G M | 3X P such that i(X Ap )uo p = (dH) p }. 

So, there exist /c-vector fields X on M which satisfy equation ([!]) on the sumbanifold C\. How- 
ever, in general these X may not be tangent to C\. Therefore, we consider the submanifold 

C 2 = {p G Ci | 3X P G T^(Ci) such that i{X Ap )u£ = (dH) p } , 

and so on. Following this process, we obtain a sequence of constraint submanifolds 

...mC 3 m...C 2 ^C 1 mM 
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where, taking into account Theorem each submanifold Cj is geometrically defined by 
Cj = {p£ C,-_i | i(Y p )(dH) p = for every Y p G (TpCj^J. 

For every j > 1, Cj is called the jth constraint submanifold. 

If we denote by £(Cj)^ the set of vector fields Y in M such that YJ, € (T p Cj)^, then one 
can obtain constraint functions defining each Cj from a local basis {Z\, . . . , Z r } of vector 
fields of X(Cj_i)^ by setting £ M = i(Z /1 )dH. 

The technical procedure to obtain these constraints is the following: 

• To obtain a local basis {Z\, . . . , Z r } of vector fields of H^=i Kerw' 4 . 

• To apply Theorem [2] to obtain a set of independent constraint functions £ M = i(Zu,)dH, 
defining C\ M. 

• To calculate X = {X\, . . . , X^), solutions to @ on C\. 

• To impose the tangency condition of X\, . . . , X^ on the constraints 

• To iterate the last item until no new constraints appear. 

This is the k-presymplectic constraint algorithm. We have two possibilities: 

• There exists an integer j > such that Cj+i = Cj = Ct. In this case, Cj is called the final 
constraint submanifold, and there exist a family of /c-vector fields X-^ = (x(, . . . ,X[) in 
M, tangent to C/, such that Q holds on Cf, that is, 

[i(x{)u; A -dH}\ Cf =0. (7) 

This is the situation which is interesting to us. 

• There exists an integer j > such that Cj = 0. This means that the equations have no 
solution on a submanifold of M. 

4 /c-symplectic field theory 
4.1 The bundle of /^-velocities 

The Lagrangian formalism of /c-symplectic field theories uses the bundle of /^-velocities of a 
manifold as phase space. First we introduce the canonical structures which this manifold is 
endowed with. 

Let T^Q = TQQ . ©TQ be the bundle of k 1 -velocities of Q, with natural projections 
t a : T\Q -» TQ and : T\Q -> Q, given in 

If (q l ) are local coordinates on !7 C Q, the induced coordinates (q l ,v A ) on (rg) _1 (L r ) are 
qr*(vig,... = q l (q) , v l A (vi q , . . . , v kq ) = v Aq {q l ) ■ 
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For Z q £ T q Q, its vertical A-lift at (v\ q , . . . , v^q) G X\Q is the vector (Z^)^- 4 , tangent to the 
fiber (tJ)- 1 ^) CT1Q, given by 



9 

then A (wi g , . . . = a* 



(«l g ,...,«fcq) 



9 dv\ 

The canonical k-tangent structure on T^Q is the set (S 1 , . . . , of tensor fields of type (1, 1) 
defined by 

S A (Wq)(Z Wq ) = (T Wq (T^(Z Wq )) V HWq) , 

for Wg = (v lq , v kq ) G T\Q, Z Wq G T Wq (TlQ). 
In coordinates we have S A = ® dg\ 

The Liouville vector field A G j£(T^Q) is the infinitesimal generator of the flow tp: R x 
T*Q — > T\Q 

ip(s; viq, ... , w fc g) = (e s t;ig, . . . , e s w fe(? ). 

Observe that A = A% + . . . + A&, where each A^ G X(T^Q) is the infinitesimal generator of 
the flow ip A : R x T^Q — ► T\Q 

1p A (s; Vlq, V kq ) = (Vl q , V( A -l)q, ^ V A q, V(A+l)q, ■ V kq) 

d 

In local coordinates we have A = Y1a=i ^ a = V% A 



dv A 



Now we want to characterize the integrable k- vector fields on T^Q such that their integral 
sections are first prolongations (jy 1 ' of maps <fi: ~R k — > Q. Remember that a k- vector field in TlQ 
is a section T : T^Q — > T^(T^Q) of the canonical projection t t iq : T^T^Q) — > T^Q. Then: 

Definition 5. ^4 second order partial differential equation (sopde) is a k-vector field T = 
(Ti, . . . ,T k ) on TfcQ which is also a section of the projection Tlr: Tu(T^Q) — > T^Q; that is, 

T\t o T = Id T i Q , 

or ; w/iai is equivalent, T Wq r ■ TA(w q ) = VAq, for w q = (vi q , . . . , V}- q ) G T^Q. 

d d 

If the local expression of the fc-vector field T = (Ta) on TlQ is T a = (X 'aYtt^ + (TaYbtt^i 

oq l ov B 

then T is a sopde iff (Ta)* = v A : 

r A (q l ,v\) = v\4j + (r A y B - d 



where {YaYb are functions locally defined in TlQ. 

If tp: R k —> T k x Q is an integral section of T = (Ti, . . . ,T k ), locally given by ip(t) = (V>*(i), V>b(*))> 
then from the last expression and Definition [3] we deduce 

~&j£t~ YAK " n dt A 
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Proposition 3. Let T = (Ti, . . . ,Tk) be an integrable SOPDE. If if) is an integral section ofT 
then tp = where <jy^ is the first prolongation of the map <f> = r o ip: R k ^ TlQ ^ Q, and 
<j) is a solution of the system of second order partial differential equations 

d 2 d/ 



Conversely, if (f) : M fc — > Q is any map satisfying (EjJ, then (j)^ is an integral section of T = 

(r 1 ,...,r k ). 

From (|HJ) we deduce that if T is an integrable SOPDE then (TaY b = (^b)a- 

Finally, using the canonical A:-tangent structure of T^Q, we have that a k- vector field T = 
(Ti, . . . , r fc ) on TlQ is a SOPDE if, and only if, S A {T A ) = A A (A fixed). 

4.2 A;-symplectic Lagrangian field theory 

Let L G C°°(T^(5) be a Lagrangian function. We define the Lagrangian forms 

e A = \s A ) o dL g q\t\q) , u£ = -&e A e &{t\q) . 

and the Lagrangian energy function 

E L = A(L) — L & C°°(Tig) . 

9L / dL \ dL 

They have local expressions 6 A = -—r- dq\ u A = dq z A d I — — ) , El = v\——. L. 

ov A \ ov A J ov A 

We introduce the Legendre map of L, which is its fibre derivative FL: T\Q — ► {T\)*Q. It 
can be defined as follows: for q G Q, u q G T q Q, (vi q , . . . , v kq ) G {T\) q Q, 

[FL(v lq , v kq )] A {u q ) = -^-L(v lq , ...,v Aq + su q , v kq )\ s=0 . 
Locally, FL(q\v A ) = ^q\ J~T^- Furthermore, we have that 6 A = FL*(9 A ), u A = FL*(uj a ). 

Definition 6. The Lagrangian L is regular if the following equivalent conditions hold: 

( d 2 L \ 

1. is everywhere nonsingular. 

\dv l A dv> B ) 

2. The second fibre derivative FL: Tj^Q — ► (T\)*Q <S> (T^)*Q is everywhere nonsingular. 

3. FL is a local diffeomorphism. 

4- {T\Q,uj a ,V = KerTrg) is a k-symplectic manifold. 



The Lagrangian L is called hyperregular if FL is a global diffeomorphism. 
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We must point out that, in field theories, the notion of regularity is not uniquely defined (for 
other approaches see, for instance, [31 17] |2"51 |2"%] l2l>] ) . 

Our purpose, however, is the study of singular Lagrangians, i.e., those which are not regular. 
Following [15], we will deal with singular Lagrangians satisfying some regularity conditions: 

Definition 7. A singular Lagrangian L is almost-regular if 

1. V := FL(T\Q) is a closed submanifold of (T\)*Q. 

2. FL is a submersion onto its image. 

3. The fibres FL~ 1 (p), for every p € V, are connected submanifolds ofT^Q. 

If L is regular, (T^Q, iv£, El) is a /c-symplectic Lagrangian system, otherwise it is a k- 
presymplectic Lagrangian system. Therefore, (T^Q, oj£, dE B ) is a fc-symplectic or a A;-presymplectic 
Hamiltonian system, depending on the regularity of L. 

In a natural chart of T^Q we have the Euler-Lagrange (EL) equations for L, which are 



d f dL 



dt A V dv\ 



dL_ 



whose solutions are maps ip: M. k — > Tj,Q that are first prolongations to Tj.Q of maps <fi = 
Tq o ip: M. k — > Q; that is, ip are holonomic. We will show that these equations can be given a 
geometric interpretation using the /c-presymplectic structure. 

Indeed, consider a map ip: M fc — > T^Q which is holonomic. Then the Euler-Lagrange equa- 
tions for ip can be also written as 

i(<p$)u£ = dE L . (10) 

As in our general discussion on fc-presymplectic Hamiltonian systems, a convenient way to 
represent the solutions of these equations can be set in terms of k- vector fields. Let us introduce 
the set 3L k L (T\Q) of A;- vector fields T = (r x , . . . , T k ) in T\Q which are solutions of 

i{Y A )u£ = dE L . (11) 

d ■ d 

If = (Y A Y-— + (r^) B — locally, then (fTTl) is equivalent to 

oq l ov l B 

{dqidv^ ~ dq3dv\) A > ~ dv\dv B B A dq^ A ~ dq i 

8 2 L _ 4 , 8 2 L 



-(T A y = — -v\. (12) 

dv B dv A dv B dv\ 

If, in addition, Y is required to be a SOPDE, i.e. (Y A y = v A , then the above equations are 
equivalent to 

d 2 L ,- d 2 L _ dL 
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These equations imply that, if T is an integrable SOPDE, its integral sections are holonomic and 
they are solutions to the EL-equations. 

If L is regular, solutions to (jlip always exist, although they are neither unique, nor necessarily 
integrable. However, if T is integrable, then the second group of equations (|12p imply that its 
integral sections are holonomic and they are solutions to the EL-equations. Hence T is a SOPDE. 

If L is not regular then, in general, equations (jlip have no solutions everywhere in T\Q but, 
in the most favourable situations, they do in a submanifold of T\,Q which is obtained by applying 
the /c-presymplectic constraint algorithm developed in Section 13.31 Nevertheless, solutions to 
equations (jll|) are not necessarily SOPDe's (unless it is required as an additional condition). In 
addition, if they are integrable, their integral sections are not necessarily holonomic, and thus 
they are not solutions to the EL-equations Q. The geometric analysis of this problem must 
be done in a separate way. (For the multisymplectic formalism of field theories, a study of this 
problem can be found in [29J). 

4.3 fc-symplectic Hamiltonian field theory 

The Hamiltonian formalism of /c-symplectic regular field theories uses the bundle of A^-covelocities 
of a manifold as phase space. 

So, consider the /c-symplectic manifold ({T\)*Q,uj a ,V), and let H € C°°((T£)*Q) be a 
Hamiltonian function. Then ((Ti,)*<2, uj a , dH) is a /c-symplectic Hamiltonian system. 

In particular, if (T\Q, oj a , d£x) is a Lagrangian system, then: 

• If L is hyper regular, we may define the Hamiltonian H = El°FL~ 1 , and ((T\)*Q, u A , dH) 
is the k-symplectic Hamiltonian system associated with L. 

• If L is almost-regular, let V be the image of the Legendre map, and j : V ^ (T£)*Q the 
corresponding embedding, and denote by FLq : T^Q — ► V the restriction of the Legendre 
map defined by jqo FLq = FL. Then, the condition of almost-regularity implies that there 
exists Hq 6 C°°(V) such that (FLq)*(Hq) = El. Furthermore, we can define uiq = Jq(uj a ). 
With these definitions, the triple (V,ui A ,dHo) is the k-presymplectic Hamiltonian system 
associated with L, and the corresponding Hamiltonian field equation @ is 

i(X° A )^ = dH 

where X° = (X®, . . . , X?) (if it exists) is a fc-vector field on V. Once again, in general, this 
equation has no solutions everywhere in V but, in the most favourable situations, they do 
in a submanifold of V which is obtained applying the /c-presymplectic constraint algorithm 
developed in Section 13.31 
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5 Applications and examples 

5.1 The Skinner— Rusk unified formalism for fc-symplectic field theory 

The so-called Skinner-Rusk formalism \A8\ [4*9] was developed in order to give a geometrical 
unified formalism for describing mechanical systems. It incorporates all the characteristics of 
Lagrangian and Hamiltonian descriptions of these systems. This formalism has been generalized 
to the fe-symplectic description of first-order field theories in [35]. Next we outline the main 
features of this formalism. 

Let us consider the direct sum TlQ © {Tl)*Q (of vector bundles over Q), with coordinates 
(q^v^pf), and denote by pry. T*Q © (T fe 1 )*Q -» T^Q and pr 2 : T^Q © (T fc 1 )*Q -> (T^)*Q the 
canonical projections. In this manifold, we have some canonical structures. 

First, if ((o»o)i, . . . , (^o)fe) is the canonical fc-symplectic structure on {T^)*Q, its pull-back 
through pr 2 yields a A;-presymplectic structure (Q±, . . . on T^Q © (T^)*Q: the 2-forms are 
defined by 9, A = (pr 2 )* (uJq) A . 

We can also define the so-called coupling function C : T^Q © (T^ )*Q — ► K by 

C(«i g , . . .,v kq ,a\, ...,a k q ):= (o^,v Aq ). 

Now, consider a Lagrangian L E C°°(r^Q). We can define a Hamiltonian function H € 
C°°{TlQ © (T fe 1 )*Q) as W = C - pr^L): 

H(v lq , . . . , v kq , a\, . . . , a q ) = C(v lq , v kq , a\, . . . ,a q ) - L(v lq , v kq ), 

which in local coordinates reads TC = afv A — L(q l ,v A ). 

Then (T^Q © (T£)*Q, £Ia, 7~i) is a A;-presymplectic Hamiltonian system; where H a=1 {Ia is 

r 8 i 

locally generated by the vector fields < — — > . We look for the solutions of its HDW equation 

[dv A ) 

which are integral sections ip: M. k — > TlQ © (T£)*Q of some integrable /c-vector field Z = 
(Z u ...,Z k )on TlQ © (TlYQ, satisfying 

iz A n A = dH. (13) 

This equation gives various kinds of information. In fact, writing locally each Za as 

Za = {Za)1 w + {ZAYB ^f B + {ZA)f ^f ' 

equation (fl3l) amounts to the following conditions: 



pt = §^°Pn, (Z A y = v A , (Z B )f = ^-o Pn . 

The first group of equations are algebraic rather than differential, and they define a sub- 
manifold Ml of T\Q ®q (T^)*Q where the equation ([13]) has solution. These constraints can 

/ 8 \ 

also be obtained by computing i — - dTC, as noted in the discussion of the fc-presymplectic 

\8v\J 
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constraint algorithm. Observe that the submanifold Ml is just the graph of the Legendre 
map FL defined by the Lagrangian L, and hence it is diffeomorphic to TlQ. We denote by 
j: M L ^TlQ 0q (T£)*Q the natural embedding. 

The second group of equations are a holonomy condition which means that the /c-velocity 
part of the integral sections of the vector field Z is the lift of a section <fi: M. k — > Q. 

The third group of equations establishes some relations among some of the coefficients (Zj^jf 
of the vector fields Za- 

Given a solution Z = (Z±, . . . , Z^) of equation (fT3|) . the vector fields Za are tangent to the 
submanifold Ml if, and only if, the functions C z i^p^ — ~~~~f °Fij vanish at the points of 

Ml, for every A,B,j (the symbol C denotes the Lie derivative). Taking into account the above 
results, this is equivalent to 

^ 3 dv B dq l A dtP B dv l c ^° ^ 

In general, equations (|13p have not a unique solution. If L is regular, taking into account 
the above results, one can define local fc-vector fields (Z\, . . . , Z^) on a neighborhood of each 
point in Ml which are solutions to (|13p . The vector field Za may be locally given by 

(ZaY = v\, {Z A )f = \^5 B A , 

with (Za)b given by equation ([H]) . Then, using a partition of the unity, one can construct global 
k- vector fields which are solutions to (|13p . When the Lagrangian L is singular one cannot assure 
the existence of consistent solutions for equation (|13p . Then, in the best cases, the constraint 
algorithm will provide a constraint submanifold Vf where these solutions exist. 



If Z is an integrable k- vector field solution to (I13p . then every integral section of Z is of 
the form ip = (i^l^h), with ipL = P r i o if;: ~R k — > T^Q, and as ip takes values in Ml then 
ipH = FL o tpL] in fact, 

Mt) = (pr 2 o VO(t) = mt), tf(t)) = ^(t), ^|^ { J = (FL o ^ L )(t) . 

Furthermore, it can be proved (see [35]) that ipL is the canonical lift of the projected section 
<fi = tq opn oil;: M. k — > Q, which is a solution to the Euler-Lagrange field equations, and that, if 
L is regular, then ipji = FL o ipL is a solution to the Hamilton-De Donder-Weyl field equations, 
where the Hamiltonian H is locally given by H o FL = El- In the almost-regular case, this last 
result also holds, but the sections ip, tpL and take values not on Ml-, F\Q and (T^)*Q, but 
in the final constraint submanifold Vf and on the projection submanifolds pri^j) T^Q and 
pr2(Vf) (T%)*Q, respectively. 

In this way, every constraint, differential equation, etc. in the unified formalism can be 

translated to the Lagrangian or the Hamiltonian formalisms by restriction to the first or the 

dL 

second factors of the product bundle. In particular, the constraint conditions pf — opn = 

ov A 

generate, by pr2-projection, the primary constraints of the Hamiltonian formalism for singular 
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Lagrangians (i.e., the image of the Legendre transformation, FL(T^,Q) C (T£)*Q), and they are 
the primary Hamiltonian constraints. 

5.2 Example 1 

Let us study a simple example, the electromagnetic field in 2 dimensions. The base manifold 
is Q = R 2 , with local coordinates {q l -,q 2 )-, and k = 2. The induced coordinates on TrjM 2 are 
(q 1 , q 2 , v }, v\, vf, v 2 ). The electromagnetic field Lagrangian is L = \{v\ + v 2 ) 2 (see [25]). 

The canonical 2-tangent structure on T3,1R 2 , (S 1 ,^ 2 ), is S 1 = ® dg 1 + <8> dg 2 and 

1 1 9 

5 2 = ® dg 1 + ® dq 2 , and the Liouville vector field reads as A = Ai + A2 = v\——. r + 
2 d 1 d 2 ^ 

Ul ^f +U2 ^ 2 r + W2 ^|- 

The Lagrangian forms are 

9\ = dL o = {v\ + vj)dq 2 , 6 2 L = dL o S 2 = {v\ + vl)dq l , 

U)\ = -del = dq 2 A dv\ + dq 2 A dvf , u| = -d0 2 L = dq 1 A dt;^ + dq 1 A du 2 , 
and the Lagrangian energy function is 

E L = A(L) — L = (v\ + v 2 ) 2 -\{v\+ v 2 ) 2 = \{v\ + v 2 ) 2 . 
Id d \ 

Since Kerw^ D Keru; 2 = ( tt~[, tt^ )> L is not regular and (T^M 2 , (w^, a; 2 ), dS^) is a 2- 



presymplectic Hamiltonian system. 
The field equation is 

i{X 1 )uj 1 L + i{X 2 )ujl = dE L 
for a 2- vector field X = (Xi,X 2 ) on T*IR 2 . 
If we write in coordinates 



_d_ 

then the field equation reads as 



(X l ) 2 (dv\+dv 2 )-{(X l )\+(X l ) 2 Mq 2 HX2)\dv^ 
We obtain that 

(X x ) 2 + (X2) 1 = v\+ v 2 , (Xx)l + (Xx) 2 = 0, (X 2 )l + (X 2 ) 2 = 0. 

Since dE^ G (Keru>^ n Keru; 2 ^) , there are no constraints and the equation has solutions at the 
whole manifold T<j,R 2 . The general solution has the form 

Xl _ ( x l}1 JL + ( „; + + B JL + C JL _ C JL + D JL, 
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where (Xi) 1 , (X 2 ) 2 , A, B, C, D, E, F and G are arbitrary functions. 
In order that the integral sections of solutions be holonomic, 

(Xi) 1 = v\, {X 2 ) 2 =v 2 2 , A = 0, C = E and D = -F, 

and, furthermore, we must also demand that [Xi,^] = 0. 

Now we will study the Hamiltonian formalism. Let (q 1 ,q 2 ,p\,p 2 ,p 2 ,p 2 ) be the induced 
coordinates on (T^)*IR 2 . The Legendre map FL: T^M 2 — ► (T^)*M 2 locally reads 

FL(q 1 ,q 2 ,v\,vl,Vi,v%) = {q\ q 2 ,p\ = 0,p 2 = v\ +v\,p\ = v\ + v\,p\ = 0). 

The image of FL, V := FL{T\Q) = {p\ = 0,p 2 = 0,pj = pi}, is a submanifold of (T^)*M 2 . 
Let Jo{q l ,q 2 ,p) = (q 1 , q 2 ,0,p,p,0) be the natural embedding and FLq: T^Q — > V the restric- 
tion of the Legendre map. We have the Hamiltonian function Hq = ^p 2 (which is such that 
(FL )*H = E L ), and the 2-forms 

ul = fooj 1 = f^dq 1 A dp\ + dq 2 A dpi) = dq 2 A dp, 

oj 2 = f co 2 = f (dq l A dp 2 + dg 2 A dp 2 ,) = dq 1 A dp. 

With these definitions, (V, oJq , lo 2 . dHo) is the 2-presymplectic Hamiltonian system associated 
with L. The corresponding Hamiltonian field equation is 

i{Y 1 )ul + i{Y 2 )J 2 = dH , 

where Y = {Y\, Y 2 ) is a 2- vector field on V. 
If, in coordinates, 



the equation is 
and we obtain 



(y) 2 dp - (Y.fdq 2 + (r 2 ) 1 dp - (Y 2 )°dq 1 = pdp 



(y 1 ) 2 + (y 2 ) 1 =p, (Y 1 )° = o, (y 2 )° = o. 

Since Ker^Q n KercjQ = {0}, there are no constraints and the equation has solutions at the 
whole manifold V . The general solution has the form 



a 

dq 2 ' 



where (Yi) 1 , (Y 2 ) 2 and A are arbitrary functions on V. 
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5.3 Example 2 

In this example we consider two independent variables (t, s) GK 2 , thus k = 2. The field compo- 
nents (dependent variables) are (q,e) E Q = W 1 x The corresponding natural coordinates 
of ffi 2 TQ are written (q,e;qt,q s ,et,e s ), and those of ® 2 T*Q are (ft e;p*,p s , 7r*, ir s ). 

We consider as Lagrangian function 

L=Y e (qt) 2 + \m 2 e-^(q s ) 2 , 

with m, t parameters, and for instance (qt) 2 is the square of qt with respect to the Euclidean 
inner product of M d . From L we compute the Lagrangian energy 

EL = Y e ^tf - \{qs? - \m 2 e 

and the Legendre map FL: © 2 TQ -> © 2 T*Q: 

FL(ft e, g s , e t , e 8 ) = (q, e, \ u -rq s , 0, 

It is clear that the primary Hamiltonian constraint submanifold Vq C © 2 T*Q is described 
by the primary hamiltonian constraints 

vr* w 0, tt s « 0. 

This also shows that the Lagrangian L is almost-regular. 
Hamiltonian formalism 

Using (q,e,p t ,p s ) as coordinates on the submanifold Vo, its 2-presymplectic structure — the 
pull-back ot the canonical 2-symplectic structure of © 2 T*Q — is given by uJq = dq A dp* and 
ujq = dq A dp* — in these expressions a summation over the invisible vector indices of q and the 
momenta is implicit. Then 

Kerc^o fl Ker^o - 



0_ 

de 



The Hamiltonian function on Vq is 



Ho = l( P t ) 2 -lm 2 e-±(p s ) 2 . 



Consider X = (X t ,X s ), a 2- vector field on Vq: 

d d d d 

Xt = Ft d~q +ft d~e +F 'dp~ t+FtS dp~ s 

Xs = Fs Wq +fs d~e +Fts dpt +FSs dp~ s 
(where the capital -F's are also vector functions). The Hamiltonian field equation for it is 
i(X t )<4 + i(X s )c4 = dH : 

F t dp l + F s dp s - {Fl + F s s )dq = ep'dp' - ^dp s + \ {{p 1 ) 2 - m 2 ) de , 
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which partly determines the coefficients of X: 

F t = ep\ F s = --p s , Ft + F* = 0, 

T 

and imposes as a consistency condition the secondary hamiltonian constraint 

\ ((p*) 2 - m 2 ) « 0; 

this can also be obtained as £ = i (jj~^J = \ ((p*) 2 — m 2 ) pa 0. 

Imposing the tangency of X to T 7 ! yields no more constraints and determines partly some 
coefficients of X: 

£ x J=p t F** i 0, C^^p'F^O, 
from which the final dynamics on V\ is given by 

* ^ dq ^ de * <9p* ' <9p s 

X - ^p s ®+f®+F t ® F f ® 
s t dq s de s dp 1 1 dp s ' 

with f t , f s , Ft arbitrary functions, and Ff, Fj arbitrary but orthogonal to p*. 

Consider the particular case of d = 1 — the q variable is just a scalar. The submanifold V\ 
is given by the constraint p l = m (or p l = —m). Then, in coordinates (q, e,p s ), the dynamics 
reads 



d , d d 

x < = me lT q + ! '3i + F 'w 

1.5 . 5 



x * = — ?V + /^ 

r de 



The analysis of the integrability of the 2- vector field X = (X t , X s ) relies on the computation of 

i A "" *J = - + ">'•) I + /• - /<) | - *r) £ ■ 

Setting it to zero determines F£ = —Tmf s and a set of two nonlinear PDEs for f s , f t : 

1 s df s df 8 

— P ~^- + fs^~ = ' 
r dq de 

df s . f dfs df s 1 s df t df t 

dq de dp s r de 

Certainly there are solutions to this equations, as for instance the one given by ft = fs = 0. 
However, it does not seem easy to give an explicit description of the whole set of these solutions. 

Finally, once one has an integrable 2-vector field X, a map ip: B 2 ^ V\, (t, s) <—> (g, e,p s ), is 
an integral section iff it satisfies 

dq _ dq _ p s de _ de _ dp s _ dp s _ 

~di = me ' ~d~s = ~~' di = }u d~s =ts, ~dt = " ~ds~ 
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Lagrangian formalism 

The Lagrangian analysis can be performed in a similar way. Let us describe it more briefly. 
Using natural coordinates (q,e;qt,et,q s ,e s ) on © 2 TQ, the 2-presymplectic structure induced 

by L is described by lo = - dq A dqt ^qt dq A de and cu s = — r dq A dq s . Then 

e e A 

t s / 9 d d d 

Kerw* fl Kerw s = — , — , e— + q t — 
\de t de t de dq t 

An arbitrary 2-vector field X = (X t ,X s ) on © 2 TQ reads 

(9 d d d d d 

Xt = F t — + j t — + F tt — h F ts — — j tt — h j ts — — , 

oq oe oqt oq s det oe s 

d d d d d d 

X s = F s — h f g — h F s f— h F ss — fst~F. I" fssT: • 
dq de dqt dq s de t de s 

If it has to satisfy the second-order condition, one has moreover 

F t = q t , F s = q s f t = e t , f s = e s . 

The field equation for it is ipQ)a/ + i{X s )oj s = dEL (maybe on a certain submanifold) . 
This determines some of the coefficients and defines just one primary Lagrangian constraint, 



X = - — 7, m . The tangency to this submanifold does not yield new constraints, and 

2 V e 2 J 

some functions in X remain arbitrary. This happens regardless of whether we impose the second- 
order condition or not, the only difference being in the number of remaining arbitrary functions. 
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